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Recently, the use of Bessel beams in evaluating the possibility of using them for a new generation 
of GPR (ground penetrating radar) systems has been considered. Therefore, an analysis of the 
propagation of Bessel beam in conducting media is worthwhile. We present here an analysis of this 
type. Specifically, for normal incidence we analyze the propagation of a Bessel beam coming from 
a perfect dielectric and impinging on a conducting medium, i.e. the propagation of a Bessel beam 
generated by refracted inhomogeneous waves. The remarkable and unexpected result is that the 
incident Bessel beam does not change its shape even when propagating in the conducting medium. 

In the last decades, electromagnetic waves have found many practical applications, including ones that are different 
from their traditional utilization in the field of communications. Among these applications, the detection and char- 
acterization of buried objects of a non-metallic nature is of particular interest. In relation to this topic, Bessel beams 
have recently been considered for evaluating the possibility of using localized waves in GPR systems. 

Bessel beams (also known as localized Bessel-X waves) are of great interest in physics because of their characteristics 
of being non diffracting beams, and due to their implications with regard to the topic of supcrluminality [T]. 

Apart from these features, their nature of localized waves may make them suitable also for other field of interest, 
As mentioned previously, the possibility of using Bessel beams in GPR apparatus has recently been considered [2J. 
This type of utilization requires a knowledge of the propagation of the beam in conducting, absorbing, media [3]. 

The aim of this note is to analyze the propagation of a Bessel beam coming from a dielectric medium (air, for 
example) and impinging at normal incidence into a conducting medium, i.e. the propagation of a Bessel beam 
generated by refracted inhomogeneous waves. We will demonstrate that, quite surprisingly, the beam does not change 
its shape when propagating in the conducting medium, but solely attenuates by going far away from the interface. 

The following analysis works in the scalar approximation, since it refers to the specific system with specific sources. 
For this system, it has been demonstrated that, among the three components of the field, there is one which is 
dominant with respect to the other two. In this situation, the scalar approximation works well, and thus its use is 
justified |4]. 

Let us begin by considering a system formed by two half-spaces, 1 and 2, as sketched in Fig. [T] 




FIG. 1: Scheme of the propagation of a plane wave coming from medium 1 and impinging on medium 2: medium 1 is a perfect 
dielectric, while medium 2 is a conducting medium. The two media (which are separated by a plane interface) are supposed to 
be homogeneous and of infinite width. 8q and 62 are the incident (real) and refracted (complex) angles, respectively. Ei, E r , 
and E t denote the incident, reflected and transmitted fields, respectively. In the dielectric medium the planes of constant 
amplitude and phase are both tilted with respect to the interface. On the contrary, in the conducting medium only the planes 
of constant phase are tilted with respect to the interface, while the planes of constant amplitude are parallel to the surface. 

A plane wave coming from medium 1 impinges on medium 2 with the incident angle 9q. Let medium 1 be a perfect 
dielectric, and medium 2 a conducting medium that is separated from medium 1 by a plane interface. The incident 
plane wave Ei in medium 1 is 
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Ei = exp [ik (ax + f3y + 72)]. (1) 

This impinging plane wave is expected to give rise in medium 2 to a plane wave E t , but of complex type. This is due 
to the fact that the refractive index n 2 in medium 2 is complex, hence an exponential decrease of the real amplitude 
is expected. Thus in medium 2 we write 



E t = Texp [ik 2 (a 2 x + (3 2 y + J2%)] = Texp [ik n2(a 2 x + (3 2 y + j 2 z)}, (2) 

where fco and k 2 (= kon 2 ) are the wave numbers in medium 1 (vacuum or air) and in medium 2, respectively; T is the 
transmission coefficient, which can be found by applying the boundary conditions at the interface between the two 
media; a, ,5,7 are the real cosine directors of the direction of propagation in medium 1, while a 2 ,/3 2 ,"/ 2 are constants 
only subject to the condition that 

a2+/3|+7i = l- (3) 

For the sake of simplicity, we have taken the amplitude of the incoming wave as being equal to one, and we have 
omitted the temporal factor exp(— iujt). In a system of spherical coordinates, the quantities a, /3,7 and a 2l f3 2 ,^f 2 may 
be written as 



a = sin 9q cos ip ol 2 — sin 9 2 cos ip 2 

j3 = sin #0 sin ip (3 2 — sin 8 2 sin i\> 2 (4) 

7 = cos 9 Q 72 = cos 6 2 , 

In principle, n 2 a 2 , n 2 (3 2 , and ip 2 could be complex quantities. However, this is not so, as can be inferred by simple 
physical considerations. 

As is well known 5, 6 , in the propagation through a plane interface into a conducting medium, the field is described 
by a system of inhomogeneous equations. Inside the conducting medium the planes of constant amplitude are parallel 
to the interface, while the planes of constant phase are tilted with respect to the same surface (see Fig. [I]). The 
direction of propagation is determined by the normal to the constant phase-planes. 

Now, for our purpose, it is crucial to note that if the planes of constant amplitude are parallel to the interface, then 
both the products 7120:2 and n 2 /3 2 in Eq. ^ have to be both real, otherwise the constant-amplitude surfaces would 
be plane, but with an inclination with respect to the interface. Moreover, since the products 7120:2 and n 2 /3 2 are real, 
tan "02 is also real, since it is the ratio of two real numbers: 



n 2 p 2 n 2 sin 9 2 sin ip 2 , E , 

tani/> 2 = = — : (5) 

n 2 02 n 2 sin 9 2 cos ip 2 

hence, ip 2 is a real angle. 

We have therefore obtained a quite unexpected result: even if the angle of refraction 2 in the conducting medium 
is complex (as expected) and the refractive index n 2 is complex, the product 712 sin# 2 is a real quantity. On this basis, 
by applying the principle of the continuity of the phase distribution at the interface (z — 0) between the dielectric 
and conducting media, we reobtain the Snell law, that is 



n 2 sin t> 2 — no sin Oq 



(6) 



from which it turns out that 9 2 is complex: sin #2 = "0 sin do /n 2 . By solving Eq. ^ we obtain the real 9 2r and 
imaginary 9 2 i parts of the complex angle of refraction 9 2 : 



sin(6» 2r ) = 
cosh(6> 2l ) = 



1 + B\n\ r + nl) + y/[l + B^nl + n^-AB^ 



2„2 
2r 



1/2 



1 + B\n\ r + n 2 2i ) + ^[l + B^{nl r + n 2 2i )Y - 4B 2 n 
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where ni r and tin are the real and imaginary parts of the refractive index ri2 , respectively, and 

B 2 = [n a sind Q /{n 2 2r + n 2 2l )] 2 . 
Coming back to Eq. ([2]), it is expedient to refer to a cylindrical coordinates system [p, <f>, z) such that 

x = pcos(f>, y — psincj), z = z. (7) 
Thus Eqs. ([I]) and ([2| can be written as (no = 1) 

Ei = exp[ifcop(sin6*o cosi/> cos <j> + sin#o siaip sin <p)\ exp(ifco cos 9qz) = 

= exp(ifco cos 6qz) exp[ifcoP s i n cos ('/ , — 4>)], (8) 

and 

E t — T exp[ifco/9(ri2 sin 62 cos ip cos <fi + n% sin 62 sin ip sin <j))\ exp(ikon2 cos 02z) = 

= Texp(ifc "-2 cos 02z) exp[ifcopsin#o cos(0 — tp)], (9) 

where use has been made of Eq. ([6]). The quantity controlling the z-dependence of the wave is complex: 



n 2 cos 6» 2 = \Jn 2 ~- sin do (10) 

which indicates that the wave attenuates in the direction normal to the interface (z-axis). 

Let us now consider a Bessel beam. As is well known, a Bessel beam originates from the interference of an infinite 
number of plane waves whose directions of propagation make the same angle, say #0 j with a given axis, say z. Let us 
suppose to have such plane waves in medium 1. Therefore, the field propagating in medium 1 is given by the sum of 
an infinite number of waves like that of Eq. ^ and results in 

E t — exp(ifc cos 9qz) / exp[ik psin8 cos(<p — tp)]dtp = 
Jo 

= 27r Jo(fcoPsin 6*o) exp(i/coz cos 6*o), (11) 

where Jo(fc psin # ) denotes the zero-order Bessel function of first kind. Similarly, by integrating Eq. Q between 
and 27T, we are able to obtain the field propagating in medium 2: 

E t = T exp(ikon 2 cos 82Z) / exp[iko p sin 8 cos((j) — ip)]dip = 

Jo 

— 27rTJo(fcopsin^o) ex p(*fco^2 cos02^) ■ (12) 



In looking at Eqs. (11) and (12), we arrive at this notable result: a Bessel beam propagating from a perfect dielectric 
medium to a conducting one does not change its shape during the passage. In other words, apart from an attenuation 
in the amplitude, no generic Bessel beam - but just the incident beam - is propagated. 

The field in the conducting medium still has the shape of the incident Bessel beam: the maxima, the minima and the 
zeros of the Bessel function do not change their positions in passing through an infinite number of conducting layers. 
The effect of the conducting media lies only in the attenuation of the amplitude which suffers some modifications: at 
z = + the amplitude is not 1 but the complex amplitude T, it then attenuates by going away from the interface (see 
Fig. [2]), due to the complex exponential factor. 

When the conductivity increases, the planes of constant phase of each wave tend to align themselves parallel to 
the interface. The effect of the conductivity is to reduce the value of the refraction angle: for large value of the 
conductivity the angle of refraction tends to zero. In this situation the Bessel beam loses its characteristic of localized 
wave and tends to become a plane wave. For infinite conductivity there is no propagation, and the field inside the 
conducting medium is zero: the electromagnetic energy is wholly dissipated by the Joule effect. 
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FIG. 2: Intensity of the incident Bessel beam at z = (a) together with the refracted Bessel beam (generated by refracted 
inhomogeneous waves) at z — 1 cm (b) and z — 2 cm (c) inside the conducting medium, as given by Eqs. (TT} and Q. The 
transmission coefficient T = 1 + (cos #o — n-2 cos 62 ) (cos So + 712 cos #2 ) ~ 1 has been found by applying the continuity conditions of 
the field and of its first derivative (with respect to the direction of propagation), to the border surface between the two media 
(2 = 0). Parameter values are: v = 9 GHz, 80 = 20°, «2 = 1.2 + i(0.12) (a reasonable value for a dry soil). 
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